Parametric Normal-Mode Splitting in Cavity Optomechanics 
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Recent experimental progress in cavity optomechanics has allowed cooling of mesoscopic mechani- 
cal oscillators via dynamic backaction provided by the parametric coupling to either an optical or an 
electrical resonator. Here we analyze the occurrence of normal-mode splitting in backaction cooling 
at high input power. We find that a hybridization of the oscillator's motion with the fluctuations 
of the driving field occurs and leads to a splitting of the mechanical and optical fluctuation spectra. 
Moreover, we find that cooling experiences a classical limitation through the cavity lifetime. 
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Introduction: Recently, cavity optomechanical systems 
that parametrically couple a driven high-frequency mode 
to a high-Q, low-frequency mechanical mode have been 
subject to increasing investigation [1 . They have been 
implemented in multiple ways. Optomechanical systems 
have been demonstrated or proposed that couple the me- 
chanical motion to an optical field directly via radiation 
pressure build up in a cavity [2, 3 s |4j|5j|6], or indirectly 
via quantum dots [7 or ions [8]. On the other hand, in 
the electromechanical domain, this has been realized or 
proposed using devices such as (superconducting) single 
electron transistors [9] ITS] . LC circuits [11], a sapphire 
parametric transducer [12], Cooper pair boxes p~3j [14], 
or a stripline microwave resonator [15 . Importantly, the 
parametric coupling can not only be used for highly sensi- 
tive readout of mechanical motion [2] but also by virtue 
of dynamical backaction be used to cool the mechani- 
cal oscillator. Indeed, recent progress has enabled the 
observation of radiation pressure dynamical backaction 
cooling [4] [5] [S] as predicted decades ago [2] ITS] . En- 
abled by this work, one emerging goal in this context is 
ground state cooling, which may open up the possibility 
of studying nonclassical states of motion or entanglement 
in mechanical objects [8j [17j [18]. For both electro- and 
optomechanical systems, it has been shown that ground 
state cooling is only possible in the resolved sideband 
regime (RSB) where the mechanical resonance frequency 
exceeds the bandwidth of the driving resonator p~9j [20] . 
This result is analogous to the laser cooling of ions in the 
"strong binding" regime [2T] . RSB cooling has recently 
been demonstrated [22[ [23] . 

Here we show that the cooling of mechanical oscilla- 
tors in the RSB regime at high driving power can entail 
the appearance of normal- mode splitting (NMS). NMS 
— the coupling of two degenerate modes with energy 
exchange taking place on a timescale faster than the de- 
coherence of each mode — is a phenomenon ubiquitous 
in both quantum and classical physics. A prominent real- 
ization occurs when atoms are coupled to a cavity field, 
which leads to the splitting of the cavity transmission 
into a doublet [24]. In addition to atom-photon inter- 



actions, NMS also arises in exciton-photon and phonon- 
photon interactions [25] . NMS has also been observed 
with "artificial atoms" in circuit QED [26] and single 
quantum dot cavity QED [27 settings. In these exam- 
ples the NMS corresponds to a splitting in the energy 
spectrum of the coupled two-mode system which may be 
accessed via linear response. In contrast, the optome- 
chanical NMS studied here involves driving two para- 
metrically coupled non-degenerate modes out of equilib- 
rium. Hence as will be discussed further below only in 
a "shifted" [19] rotating-frame representation does the 
Hamiltonian become analogous to the one characterizing 
the aforementioned examples. Concomitantly, the split- 
ting, rather than appearing directly in the cavity trans- 
mission, manifests itself in the fluctuation spectra. This 
scenario is reminiscent of the single trapped ion realiza- 
tion of the Jaynes-Cummings model [28] with the role of 
the pseudospin now played by the optical (or electrical) 
mode. Since this type of normal-mode splitting occurs 
during RSB cooling, we analyze how the onset of NMS 
affects and limits cooling in the RSB regime. 

Theoretical model: We start from the rotating- 
frame Hamiltonian H' = —h/^'a^a p + ^ m aJ n a m + 
^^ m aja p (a m + a^) + h(s+a p + s+aj) which provides 
a unified treatment of both a coherently driven opti- 
cal and electrical resonator (frequency uo p ) coupled to 
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Figure 1: (a) Electromechanical realization of parametric cou- 
pling of a mechanical oscillator to an LC circuit, where the 
coupling is determined by ^fU=o = 2dc t lt \ x=0 ^ tot is the 
total capacitance), (b) Optomechanical realization of para- 
metric coupling of a mechanical oscillator to a Fabry- Perot op- 
tical mode with ^^-\ x=0 = — ^-\ x =o (L is the cavity length). 
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a mechanical oscillator (frequency Q m <C co p ) via the 
dimensionless parameter 77 = (xo/Q m )-^-\ x= o. Here 
xq = ^h/2m e fffl rn is the zero point motion of the me- 
chanical mode, m e ff its effective mass, A' the detuning 
of the drive from co pi and a m (a p ) is the annihilation oper- 
ator for the mechanical (optical or electrical) mode. The 
dependence of the resonant frequency co p on the mechani- 
cal oscillator's deflection x determines the strength of the 
coupling via ^-\ x=0 [cf. Fig.jl]. The driving rate is given 
by |s_|_| = ^P/Hco p r exj where P denotes the launched in- 
put power and r~ x x is the external coupling rate. 
We derive the Heisenberg equations of motion for the 
canonical variables and introduce noise operators £ m (t) 
and £ p (t) weighted with the rates T m and ft that charac- 
terize, respectively, the dissipation of the mechanical and 
optical (or electrical) degree of freedom. Subsequently, 
we shift the canonical variables to their steady-state val- 
ues (i.e. a p — > a + a p and a m — > /3 + a m ) and linearize 
to obtain the following Heisenberg-Langevin equations 

HHIEQIES]: 

a p = (iA -^) a p- { a m + a m ) + y/K£ p (t) , (1) 

Here, A is the detuning with respect to the renormalized 
resonance and A < leads to cooling [19]. The optome- 
chanical coupling rate is given by g m = 2ar]Vt rni which 
is positive by an appropriate choice for the phase of s + , 
and \a\ 2 gives the mean resonator occupation number. 
In the case of the mechanical degree of freedom, the 
rotating wave approximation in the coupling to its envi- 
ronment implied by Eqs. ([!]) is only warranted for high 
Q values (and small g m /^m) HE] — conditions that are 
satisfied in the parameter regime of interest for ground 
state cooling. The latter also requires T m <C ft, which 
we will assume throughout our treatment. Equations ([I]) 
and their Hermitian conjugates constitute a system of 
four first-order coupled operator equations, for which 
the Routh-Hurwitz criterion implies that the system is 
only stable for g m < ^(A 2 + ft 2 /4) Q m /\A\ w ^ m (if 
Vt rn ^> ft and |A| w f2 m ). 

Here, we follow a semi-classical theory by considering 
noncommuting noise operators for the input field, i.e., 

<&(*)> = o, <£(*%(*)> = *hW - *), (Ut'Kit)) = 

(n p + l)5(t f — t), and a classical thermal noise in- 
put for the mechanical oscillator, i.e. = 0, 
<&(f )&»(*)> = (UV)tL(t)) = n m 5(t' - t), in Eqs. §. 
The quantities n m and n p are the equilibrium occupation 
numbers for the mechanical and optical (or electrical) os- 
cillators, respectively. We transform to the quadratures 
(i.e., x/xq = a m + aj^) and solve the Langevin equations 
in Fourier space [20 j. Thus we recover a steady-state 
displacement spectrum [29] given (for n p = 0) by 



x -\u) = n 2 m + 2n m n s (uo) -co 2 -ico [r m + v s (co)\ (2) 
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This spectrum is characterized by a mechanical suscep- 
tibility x(^) that is driven by thermal noise (ex n m ) and 
by the quantum fluctuations of the radiation pressure 
(quantum backaction). In linear cooling theory the 
susceptibility is approximated by evaluating the terms 
T s (co) and Q s (lj) at the (bare) mechanical frequency 
[4j [6] [30]. Then T s (^ m ) coincides with the cooling rate 
and is linear in the input power (g^ cx P). 

Parametric NMS: The above approximation is only 
adequate for weak driving such that g m <C ft [19, 20]. To 
obtain an understanding of the mechanical susceptibility 
beyond this linear regime, we return to the linearized 
Heisenberg-Langevin equations ([I]) and calculate the 
corresponding eigenfrequencies that determine the dy- 
namics of the system. Though there exists an analytical 
solution, it is rather opaque and does not provide physi- 
cal insight, so that we will use instead an approximation 
scheme appropriate for the parameter regime relevant 
for the observation of NMS and to attain ground state 
cooling. Along these lines, we focus in the following on: 
(i) the RSB regime (ft < Q m /2) necessary for ground 
state cooling [191 [20j [22], (ii) optomechanical coupling 
9m % ^m/2, and (in) 5 2 <C Q 2 n (S = —A — Q m , the fre- 
quency detuning from the lower sideband) . In the shifted 
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Figure 2: Real and imaginary parts of the eigenvalues [cf. 
Eqs. (J3| and Q] of the linearized cooling problem corre- 
sponding, respectively, to the eigenfrequency and mode damp- 
ing for A = — Q m and ft/O m = 0.2. The inset magnifies the 
resonance shift before the mode splitting. The real part is un- 
derlaid with the normalized classical displacement spectrum 
(contribution oc n m ) [cf. Eq. [5], thereof sample curves are 
highlighted in Fig. [3] for the rates marked by the dots. 
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representation corresponding to Eqs. ([!]) [19] the relevant 
part of the parametric interaction in Hamiltonian H' 
is described by an effective dipole-like interaction term, 



i.e. hr)Q m aj ) a p (a ri 



) _> *f±(a p + al)(a rt 



after neglecting the nonlinear term. This interaction 
term is analogous to the Jaynes-Cummings setting 
(with a p — > a - ) and naturally leads to resonance 
splitting when the modes have matching frequencies. 
The off-resonant counterrotating terms (CRT) ex aja^, 
a p a m induce a small frequency shift analogous to the 
Bloch-Siegert shift in atomic physics [31 . These CRT 
terms, which are responsible for the mixing between the 
creation and annihilation operators in the Heisenberg- 
Langevin equations ([!]), can be treated in perturbation 
theory within the parameter range defined by (i)-(m). 
The first nonvanishing order in this perturbative expan- 
sion is quadratic in the CRT and yields a correction to 



the decoupled eigenvalues u± 

(2) 

we take T m = in uj±]: 



,(°) 



,( 2 ) 



[note that 



s ft + rv, 

2~ l — 



,(2) 



±^V<&-(«/2-r m /2 + i<5)2, 

9]J± 

2tl m + 5± ytgl- ( K /2 + iSy) 



(3) 
(4) 



Naturally, there is another pair of eigenfrequencies given 
by —u J. In Fig. |2| the real and imaginary parts of the 
eigenvalues are plotted. The inset shows the frequency 

(2) 

shift due to the CRT. If we choose the value S = 

(i.e., A = — Q m ) relevant for n <C ^ m (see below) and 
neglect T m , the square root term of leads to two 
regimes. While for g m < k/2 the term is fully imaginary 
and modifies the decay rate of the modes, for g m > n/2 
it becomes real instead and the real parts of the eigenfre- 
quencies exhibit the splitting that signals NMS (Fig. [3|. 
The latter is associated to a mixing between the mechan- 
ical mode and the fluctuation around the steady-state of 
the resonator field. Classically, this fluctuation can be 
understood as a beat of the pump photons with the pho- 
tons scattered on resonance which leads to oscillations 
with frequency |A| in the intensity time-averaged over 
2tt/co p . For n 2 /A^g^ the splitting (~g m ) is proportional 
to the square root of the mean cavity photon number 
(a 2 ). This is analogous to NMS in atomic physics where 
the splitting of the cavity resonance is proportional to the 
square root of the number of atoms coupled to the cavity 
mode [24 . When detecting the phase fluctuations in the 
transmitted light with a homodyne detection scheme, the 
signal at Q m splits [cf. Fig.^], but the (suppressed) scat- 
tered light at the carrier frequency exhibits no splitting. 
It is important to note that the splitting in the displace- 
ment spectrum is not observed unless g m > Kj\[2 due to 
the finite width of the peaks. Due to the requirements 



on the cavity bandwidth and the detuning, the param- 
eter regime in which NMS may appear implies cooling. 
In turn, for a positive detuning (which entails amplifica- 
tion) the observation of NMS is prevented by the onset 
of the parametric instability [3]. Therefore, a discussion 
of NMS cannot be decoupled from an analysis of the as- 
sociated cooling. We also show below that the CRT in 
the interaction lead to the quantum limit of backaction 
cooling QH |U. 

Effect of NMS on backaction cooling: We now use 
the approximate eigenfrequencies to perform contour in- 
tegration on the normal ordered mechanical spectrum in 
order to obtain the final occupancy of the mechanical os- 
cillator rif = (a] n (r)a m (0)) | - In this treatment we 
take both the thermal and the vacuum noise of the driv- 
ing resonator into account. A finite value for n p may be 
relevant for electromechanical systems [Til IE] • Within 
our approximation scheme we can introduce a formal pa- 
rameter that tags the CRT terms and expand rif in its 
powers. To zeroth order the poles are determined by 
the approximate eigenfrequencies qj±\—uj±^* given in 

Eqs. (JiJ), and it is straightforward to evaluate (in- 
cluding r m ). To second order we use instead the poles 
+uj± \ —lu±^*—lu±^*. Subsequently, is expanded 
in the small parameters g m /^ m , ft/^ m , and \5\/Vt m up 
to second order with T m — > 0. Both rSp and n^p do not 
contain terms linear in (5, allowing one to directly mini- 
mize the result with respect to S by setting S — > 0. This 
yields 



,(°) 



,( 2 ) 



« 4^2 



The final occupancy rif 

contributions. One is proportional to the occupancy 



•r 



1 

2 

,(2) 



J2 



(5) 



rif ' consists of three 



mm 1m 



0.5 0.75 1.0 1.25 1.5 
Normalized Frequency C0/£2 m 



! O.Ol 



2 



b i 


k 





0.6 0.8 1.0 1.2 1.4 1.6 
Normalized Frequency C0/Q m 



Figure 3: (a) Normalized logarithm of the classical displace- 
ment spectrum (contribution oc n m ) [cf. Eq.[2] as a function of 
the normalized detuning for g m /^m = 0.4 and n/Q>m — 0.2. 
(b) Typical displacement spectra for coupling rates that are 
represented by the dots in Fig. [2] The solid curves above 
correspond to the phase spectral density S&(uj) measured in 
homodyne detection. 
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of the thermal bath n m and displays linear cooling for 
r m < g m < ^ i.e., n/ « ^fe n ™- When # m approaches 
ft, deviations from the linear cooling regime become ap- 
parent. Indeed, the final occupancy is always limited by 



n f 



> 



which implies that the largest temperature 



reduction is bound by the cavity decay rate ft [33] . This is 
equivalent to the condition Q m > n m ^ for ground state 
cooling. It is noted that operation in the deeply RSB 
regime is advantageous to avoid photon-induced heating 
[22] . entailing that the condition on the mechanical Q is 
therefore more stringent. A second contribution is pro- 
portional to the finite occupancy of the driving circuit 
(n p ) and corresponds to heating from thermal noise in 
its input. It implies that it is impossible to cool below 
the equilibrium occupation of the resonator. If we assume 
that the mechanical and electromagnetic baths are at the 
same temperature T m , it entails rif > n m ^ 1L . Last, there 

(2) 

is a term in ' that is temperature-independent and 
corresponds to heating from quantum backaction noise. 
This term determines the quantum limit to the final oc- 
cupancy and agrees with Refs. [I9j [20]. Interestingly, in 
the present analysis the quantum limit arises from the 
CRT. We note that the trade-off between the quantum 
limit and the cavity bandwidth limitation leads to an op- 
timal value for ft. Consistent results are obtained with a 
covariance matrix approach [32] . 

Finally, we consider appreciable cooling [rif <C n m so 
that we can take T m — > in the denominator of Eqs. (p5|)] 



and optimize n 



(0) 

/ 



(2) 



with respect to g mi which yields 




2 r m K(n p + 1/2) 



(6) 



for # opt = ^/4n m r m ft^/[n p + 1 /2 + n m r m /ft]). In the 
ground state cooling regime, the first three terms of 
Eq. ([6| always give the correct order of magnitude. Thus, 
a comparison of g opt with the condition g m > ft/2 im- 
plies that optimal ground state cooling leads to NMS 
only when the thermal noise [first term in Eq. (|6|] is at 
least comparable to the quantum backaction noise [third 
term in Eq. This is likely to be the case in current 
endeavors to reach the ground state. 

Experimental realization: To demonstrate that the 
observation of parametric NMS is within experimen- 
tal reach, we discuss the parameters from Ref. [22] : 
n m /27r = 73.5 MHz, k/2tt = 3.2 MHz, r m /27r = 1.3 kHz. 
The cooling rate r c /27r = 1.56 MHz was extracted from 
the displacement spectrum's FWHM. A comparison with 
Eqs. ([2| then yields a coupling rate g m /27r ~ 2.0 MHz. 
Therefore, the observation of parametric NMS is within 
experimental reach. In the electromechanical domain, 
using a superconducting coplanar waveguide resonator, 
Ref. [23] reports coupling rates of g m /27T = 6 kHz for a 
cavity with a decay rate k/2tt = 230 kHz. 
In summary, we have analyzed a novel instance of NMS 



that occurs in cavity optomechanics due to the coupling 
between the fluctuations of the cavity field and the me- 
chanical oscillator mode. Furthermore, we have eluci- 
dated its implications for ground state cooling, namely, 
the limitation through the cavity bandwidth. 
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